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Introduction 

Silicon  carbide  is  the  only  compound  species  that  exists  in  the  solid 
state  in  the  Si -C  system  and  can  occur  in  the  cubic  (C),  hexagonal  (H)  or 
rhombohedral  (R)  structures.  It  is  also  classified  as  existing  in  the 
beta  and  alpha  modifications.  The  beta,  or  cubic,  form  crystallizes  in 
the  zincblende  or  sphalerite  structure;  whereas,  a  large  number  (approxi¬ 
mately  140)  of  the  alpha  occur  in  the  hexagonal  or  rhombohedral  forms  known 
as  polytypes.  At  one  time  it  was  thought  that  the  beta  modification  was 
thermodynamically  stable  at  low  temperatures  and  that  it  transformed 
irreversibly  to  one  or  several  of  the  alpha  forms  by  a  solid  state  phase 

transformation  at  about  2100°cJ  Above  this  temperature  the  alpha  modifi- 

2 

cation  was  considered  stable.  However,  later  efforts  by  Ryan,  et  al.  and 
other  international  groups  involved  in  thin  film  growth  have  shown  that  a 
relatively  moderate  excess  pressure  (~  10  atm.),  especially  of  N2»  favors 
the  growth  of  3  at  all  temperatures.  Indeed,  under  certain  conditions, 
the  transformation  is  completely  reversible.  Impurities,  defects,  temper¬ 
ature,  Si  evaporation  and  ambient  gas  appear  to  influence  not  only  the 
transformation  but  the  resulting  polytype  as  well.  This,  coupled  with 
small  differences  in  stacking  fault  energies  among  the  various  forms, 
makes  reproducible  growth  of  thin  films  difficult  unless  considerable  pre¬ 
cautions  are  taken. 

Because  of  the  emerging  need  for  high  temperature,  high  frequency  and 

high  power  electric  devices,  blue  L.E.D.s,  Schottky  diodes,  U.V.  radiation 

detectors,  high  temperature  photocells  and  heteroj unction  devices,  silicon 

carbide  is  being  increasingly  examined  throughout  the  world  for  employment 

as  a  candidate  material  in  these  specialized  applications.  The  electron 

Hall  mobility  of  high  purity  undoped  3-Si C  is  approximately  a  factor  of  three 
2 

larger  (~  1000  cm  /v-sec)  than  that  of  the  a- form  over  the  temperature 
range  of  300-1000K  because  of  the  smaller  amount  of  phonon  scattering  in 
the  cubic  material.  The  energy  gap  is  also  less  in  the  6  form  (2.3  eV) 
compared  to  the  a-forms  (e.g.,  6H  =  2.86  eV).  Thus,  the  8-form  is  now 
considered  more  desirable  for  electronic  device  applications,  and,  therefore, 
the  growth  of  thin  films  of  this  material  constitutes  the  principle 
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objective  of  this  research  program.  Unfortunately,  the  earlier  push  in  the 
1956-1970  time  span  to  develop  SiC  as  an  electronic  material  concentrated 
heavily  on  high  temperature  growth  processes  such  as  the  Lely  sublimation- 
condensation  technique  which  produced  a  variety  of  a  polytypes  in  experi¬ 
ments  which  were  rarely  reproducible.  Toward  the  end  of  this  initial 
thrust,  techniques  such  as  chemical  vapor  deposition  (CVD),  sputtering, 
traveling  solvent  and  solution  growth  showed  promise  not  only  as  techniques 
per  se  but  as  experimental  avenues  wherein  the  growth  of  tf-SiC  could  be 
achieved. 

The  objective  of  the  present  research  program  is  to  capitalize  on  and 
extend  the  knowledge  of  the  CVD  and  sputtering  processes  that  were  discerned 
not  only  in  the  brief  initial  efforts  with  SiC,  but  that  have  also  been 
developed  for  these  processes  in  the  intervening  years  in  research  on  other 
semiconductor  materials.  Although  the  actual  experimental  portion  of  the 
effort  is  still  awaiting  delivery  of  several  components,*  considerable 
progress  has  been  made  in  terms  of  design  of  the  CVD  equipment  and  the 
computer-aided  production  of  CVD  phase  diagrams.  A  sputtering  unit  has 
arrived  just  prior  to  this  writing  and  the  peripheral  equipment  for  this 
unit  is  being  installed  as  it  is  received.  The  details  of  these  efforts 
are  reported  below. 


II.  Chemical  Vapor  Deposition 
A.  Background 

Chemical  vapor  deposition  is  a  material  synthesis  method  in  which  the 
constitutents  of  the  vapor  phase  react  to  form  a  solid  film  at  some  surface. 
Typical  gas  combinations  which  have  been  employed  in  the  growth  of  SiC  are 
CCl^  and  Si  Cl 4 ,  Si and  C^Hg,  SiCl^  and  CgH-^  and  ChgSiClg  alone.  In  private 
comnuni cations,  Ryan-*  has  informed  the  principal  investigator  that  a  small 
amount  of  HC1  is  desired  to  etch  away  the  areas  of  poor  growth  which  begin 
on  the  substrate. 


it 

The  recent  tremendous  surge  in  semiconductor  manufacturing  by  all  common 
techniques  has  generated  order  backlogs  for  equipment  and  parts  having 
waiting  periods  of  up  to  six  months,  even  for  most  favored  customers. 
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The  reaction  in  the  vicinity  of  the  heated  substrate  has  been  shown  to 

4  5 

vary  according  to  the  type  of  source  species  and  ambient  gas.  ’  Gatos 
and  co-workers^  have  analyzed  the  equilibrium  partial  pressures  in  the 
hydrogen-silicon-carbon  system  and  concluded  that  the  efficiency  of  SiC 
formation  decreases  with  1}  increasing  temperature  1650°C)  which  causes 
decomposition  and  loss  of  Si  and  2)  a  decreasing  sil icon-to-carbon  ratio  in 
the  gas  mixture.  The  epitaxial  growth  from  a  silane  and  propane  mixture 
was  found  by  these  researchers^  to  be  linearly  dependent  on  the  silane 
pressure  for  a  given  temperature  and  propane  input  which  suggests  that  Si 
mobility  and  concentration  at  the  substrate  controls  growth.  As  noted  above, 
Mogab  and  Leamy7  found  similar  results  in  their  studies  of  the  reactions  of 
hydrocarbons  with  silicon  single  crystal  surfaces.  The  kinetics  research  of 
the  present  writer  has  also  shown  that  Si  is  a  much  slower  diffusing  species 
than  C  and  should  control  growth  where  mobility  across  surfaces  is  of 
importance.  Finally  Gatos  et  al.®  have  reported  that  the  SiC  growth  process 
is  best  described  by  an  adsorption-desorption  process  and  must  be  carried 
out  in  an  atmosphere  having  a  Si/C  ratio  of  1/1  or  larger  if  the  simul¬ 
taneous  deposition  of  graphite  and  SiC  is  to  be  avoided.  These  findings 

4  8 

are  similar  to  those  of  Brander.  ’ 

In  order  to  understand  CVD  processes,  one  should  know  which  chemical 
reactions  occur  in  the  reactor  and  to  what  extent.  Furthermore,  the  effects 
of  process  variables  such  as  temperature,  input  concentrations  and  flow  rates 
on  these  reactions  must  be  understood.  The  approaches  taken  to  date  to 
discern  the  aforementioned  variables  have  been  1)  visits  and  discussions 
with  various  former  and  present  investigators  in  the  CVD  field,  2)  design  of 
the  experimental  apparatus  to  be  used  in  the  present  work  and  3)  development 
via  computer  of  CVD  phase  diagrams.  Visits  to  C.  Ryan  and  co-workers  at 
Rome  Air  Defense  Conmand,  Berkman  and  co-workers  at  RCA,  several  laboratories 
in  California  and  the  International  CVD  Conference  in  October  of  1979  have 
been  invaluable  in  the  development  of  the  design  of  the  CVD  process,  the 
determination  of  the  pertinent  equipment  for  each  major  section  and  the 
handling  of  toxic  and  flammable  gases.  Also,  these  visits  have  reinforced 
the  ideas  espoused  much  earlier  by  workers  in  the  SiC  growth  field  concerning 
the  importance  of  careful  control  of  temperature,  gas  composition  and  flow 
rate  and  the  maintenance  of  very  low  levels  of  impurities. 


ft 
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B.  CVD  Phase  Diagram 

As  the  design  of  the  CVD  system  evolved,  it  became  obvious  that  more 
information  was  needed  in  terms  of  the  solid  reaction  products  which  would 
occur  under  prescribed  conditions  of  concentrations  of  certain  reactant 
gases,  pressure  and  temperature.  Our  research  has  revealed  that  numerous 
methods  have  been  reported  for  computing  thermodynamic  equilibrium  compo¬ 
sitions  and  that  essentially  all  the  techniques  have  been  based  on  the 
methods  of  either  Brinkley  (9,  10)  or  White,  et  al.  (11).  The  former 
requires  the  writing  of  a  series  of  chemical  reactions  and  the  computation 
of  their  respective  equilibrium  constants;  whereas,  the  latter  involves 
minimizing  the  summation  of  the  free  energies  of  all  species.  Eriksson  (12) 
has  extended  White  et  al.'s  method  to  include  systems  containing  more  than 
one  condensed  phase  and  has  developed  a  computer  program,  SOL GAS,  for  per¬ 
forming  the  calculations.  He  has  modified  the  program  (13,  14)  to  be 
SOLGASMIX  so  that  condensed  solution  phases  could  also  be  considered. 

To  prevent  the  omission  of  an  important  species,  the  programs  can  consider 
all  conceivable  gaseous  species  and  condensed  phases  of  invariant  or  variable 
stoichiometry  and  mixtures  at  constant  total  pressure  and  temperature  in  a 
single  calculation.  Unstable  condensed  phases  are  rejected  from  the  calcu¬ 
lation.  This  program  can  also  handle  non-ideal  phases,  provided  activity 
coefficient  relationships  are  available.  Finally,  Bessman  and  Spear15a’  b 
have  produced  a  new  version  of  the  program  (SOLGASMIX-PV)  incorporating  the 
ideal  gas  law.  This  last  version  is  now  capable  of  calculating  equilibra 
at  a  constant  total  gas  volume  with  variable  total  pressure.  These  investi¬ 
gators  have  kindly  provided  these  programs  and  helpful  suggestions  for  use 
in  developing  CVD  phase  diagrams  in  this  Si C  research.  A  listing  of  the 
programs  as  used  in  our  research  are  given  in  Appendix  I. 

In  developing  CVD  phase  diagrams,  thermodynamic  calculations  using  CVD 
variables  can  be  used  to  show  the  ranges  of  input  conditions  for  a  CVD 
system  which  will  produce  a  specific  condensed  phase(s)  and  phase  boundaries 
at  equilibrium.  In  our  particular  employment  of  these  CVD  programs,  all 
conceivable  gaseous  species  are  included  with  the  condensed  phases  In  calcu¬ 
lating  the  equilibrium  composition  of  the  system  as  a  systematic  function  of 
the  CVD  variables  of  temperature,  total  pressure  and  reactant  concentrations. 
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Both  a-SiC  and  3-Si C  are  included  in  our  calculations.  In  each  calculation, 
one  of  the  variables  is  changed  while  all  other  variables  are  held  constant 
at  an  arbitrary  set  of  "standard  values."  These  standard  values  are  typi¬ 
cally  chosen  near  the  mid-range  of  the  experimental  controls,  so  the  desired 
condensed  phase  is  deposited  when  all  variables  are  fixed  at  their  respective 
standard  value. 

In  our  present  calculations,  the  elements  Si,  C  and  ^  are  the  principal 
species,  since  Si and  CH^  are  the  reactant  gases  and  ^  is  the  carrier  gas. 
The  type  of  substrate  becomes  important  only  if  it  emits  a  gaseous  species 
or  enters  the  reaction  upon  the  formation  of  the  condensed  phase  or  causes 
structural  rearrangement  of  the  solid  products.  Although  this  aspect  of  the 
research  is  in  its  initial  stages,  one  of  the  important  diagrams  developed 
thus  far  for  the  Si-C-H  system  at  one  atmosphere  pressure  is  presented  in 
Figure  1.  An  inert  substrate  is  assumed.  The  Si/Si+C  mole  ratio  depends  on 
not  only  the  fixed  input  gas  ratio,  but  also  on  how  much  Si C  is  deposited. 

The  same  is  true  for  the  amount  of  Hg  in  the  system.  Therefore,  a  maximum  of 
three  phases  (a  gas  phase  and  two  condensed  phases)  can  exist  in  equilibrium 
at  any  given  temperature. 

The  resulting  diagram  of  Figure  1  reveals  that  above  1000K  (the  lowest 
temperature  yet  examined,  B-SiC  has  a  negative  standard  free  energy  of  for¬ 
mation  and  the  amount  that  forms  as  a  single  phase  to  the  left  of  the  Si /Si 
+  C  =  0.5  point  occurs  over  a  varying  range  of  these  ratios  which  is  a 
direct  function  of  the  amount  of  present.  At  all  other  ratios  or  temper¬ 
atures,  solid  free  C  or  Si  also  deposit  or  Si  becomes  a  liquid.  It  should 
be  kept  in  mind  that  a  positive  AG^  for  a-SiC  does  not  mean  that  it  will  not 
form  (especially  if  a-SiC  single  crystals  were  used  as  a  substrate)  only 
that  it  is  slightly  less  thermodynamically  preferable  than  6-Si C  under  the 
given  conditions.  Since  the  thermal  decomposition  of  Si is  essentially 
complete  at  temperatures  greater  than  1000°K,  the  effect  of  adding  ^  to 
the  system  is  to  promote  the  reaction  with  excess  carbon  atoms  to  produce 
various  gaseous  hydrocarbon  species  and  thus  provide  for  the  formation  of 
6-Si C  over  a  wide  range  of  Si/Si+C  ratios  less  than  0.5. 


u-H  System  at  one  Atmosphere  Pressure 
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It  has  been  reported  by  Ban^  that  more  than  80‘,y  of  the  SiH^  in  a  2: 

Si mixture  in  H2  is  decomposed  at  T  >  1073K.  Thus,  at  lower  concen¬ 
trations  of  even  greater  decomposition  is  to  be  expected,  Based  on 
our  calculations,  almost  all  Si should  be  decomposed  under  equilibrium 
conditions.  Furthermore,  the  higher  the  temperature,  the  greater  will  be 
the  decomposition  of  CH^.  As  noted  above,  the  pressure  of  H2  will  react 
with  the  C  and  thus  provide  considerable  process  control. 

Finally  it  should  be  noted  that  diagrams  were  also  derived  using  Ar 
at  one  atmosphere  in  lieu  of  H2  from  the  foregoing  discussion,  the  phase 
field  of  3-SiC  is  essentially  that  labeled  (1)  in  Figure  1,  as  there  is  no 
species  to  react  with  the  free  carbon. 

The  next  calculations  will  entail  going  to  reduced  pressures  in  order 
to  determine  any  changes  in  amount  of  6-Si C  that  will  form.  Preliminary 
research  at  10  atm  without  the  presence  of  ^  indicates  the  formation  of 
6-Si C  only  at  a  ratio  of  Si/Si+C  of  0.5  until  approximately  1550K  where  the 
region  widens  rapidly  toward  the  lower  ratios  as  temperature  is  increased. 
These  investigations  will  continue  throughout  the  next  year. 

C.  Design  of  the  CVD  System 

At  normal  atmospheric  pressure,  the  rates  of  mass  transfer  to  the 
gaseous  reactant  and  by  product  species  and  the  surface  reactant  rate  to 
form  the  film  deposit  are  generally  of  the  same  order  of  magnitude. 

Therefore,  in  order  to  attain  the  objective  of  uniform  film  thickness  and 
properties  over  a  surface,  it  is  necessary  to  ensure  that  careful  consider¬ 
ation  be  given  to  both  types  of  rate  determining  steps.  Mass  transfer 
rates  depend  mainly  upon  reactant  concentration,  diffusivity  and  boundary 
layer  thickness  which  is  related  to  reactor  configuration,  flow  velocity, 
distances  from  edges,  etc.  Surface  reaction  rates,  however,  depend  mainly 
upon  reactant  concentration  and  temperature. 

The  diffusivity  of  a  gas  is  inversely  related  to  pressure.  Therefore, 
as  pressure  inside  the  reactor  is  lowered  from  one  atmosphere  to  0.5-1 .0  torr, 
the  diffusivity  increases  by  a  factor  of  1000.  This  is  only  partially  off¬ 
set  by  the  fact  that  the  laminar  boundary  layer  (distance  across  which  the  re¬ 
actants  must  diffuse)  increases  at  less  than  the  square  root  of  the  pressure. 
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The  net  effect  is  more  than  one  order  of  magnitude  increase  of  the  gas 
phase  transfer  of  reactants  to  and  by-products  from  the  substrate  surfac'*, 
and  the  rate-determining  step  in  the  sequence  is  the  surface  reaction. 
Therefore,  in  such  a  low-pressure  operation  (LPCVD),  much  less  attention 
need  be  paid  to  the  mass  transfer  variables  that  are  so  critical  at 
atmospheric  pressure. 

An  additional  advantage  of  LPCVD  is  the  considerable  reduction  in 
autodoping  from  gaseous  species  derived  from  the  substrate  material.  In  CVD 
processes  at  one  atmosphere,  the  laminar  portion  of  the  boundary  layer  next 
to  the  sample  is  very  narrow.  As  such,  any  species  diffusing  out  of  the 
substrate  enters  this  narrow  layer  and  has  a  relatively  high  probability  of 
returning  to  surface  of  the  material,  being  deposited  and  becoming  incorpo¬ 
rated  into  this  deposit.  For  these  reasons,  an  LPCVD  system  has  been 
designed  as  shown  in  Figures  2-5  and  described  below.  Because  this  system 
is  and  will  be  continually  in  the  process  of  evolution  even  after  the 
completion  of  the  initial  "best"  arrangement,  the  author  has  chosen  to 
present  the  figures  in  an  unfinished  form  rather  than  to  take  several  days 
to  make  engineering  drawings.  The  latter  will  be  produced  for  that  ar¬ 
rangement  which  will  be  found  to  be  optimum. 

As  shown  in  Figure  2,  the  principal  reactant  gases  are  methane  (CH^) 
and  silane  (Si )  with  a  carrier  gas  of  H2.  Although  all  gases  will  be  of 
highest  available  purity  (H2  will  be  purified  "in-house"  using  a  Pd/Ag 
cell),  CH4  and  H2  will  also  be  "scrubed"  to  remove  any  residual  oxygen  and 
H20  which  could  react  with  Si to  form  Si02>  In  addition,  crosspurge 
assemblies  will  be  employed  to  eliminate  impurity  capture  during  tank 
exchange.  HC1  will  be  used  to  clean  the  single  crystal  Si  substrates  and 
to  etch  away  any  excess  Si  from  the  SiC  during  growth.  Stainless  steel 
tubing,  specially  cleaned  for  chromatographic  experiments,  will  be  used 
throughout  the  system.  An  MKS  electronic  system  will  be  used  to  control 
the  flow  of  each  of  the  gases  in  the  system.  The  flow  of  the  carrier  gas 
(when  used)  will  be  set  at  a  particular  value  and  Si ahd  CH^  slaved  to 
this  value  in  a  particular  ratio.  HCl  will  have  an  independent  flow  con¬ 
troller,  since  1)  the  operator  may  desire  to  vary  the  amount  of  etch  used 
for  the  Si  and  SiC,  respectively,  and  2)  the  MKS  unit  can  handle  only  three 
gases.  The  flow  rate  of  each  gas  (except  HCl)  will  be  measured  by  a  Tylan 
flow  meter  which  returns  a  signal  to  the  MKS  unit  which,  in  turn,  compares 


Schematic  of  Gas  Flow,  Monitoring  and  Regulation  for  CVD  Systei 


Figure  5. 


Schematic  of 
Pumping  Systems 
for  CVD  Process 
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this  signal  with  the  electrical  setpoints  established  by  the  operator  and 
adjusts  the  Brooks  electronic  valves,  if  the  desired  and  the  actual  flow 
rates  are  different.  Check  valves  will  be  employed  in  each  line  to  prevent 
contamination  of  the  lines  and  gas  bottles.  In  addition,  argon-operated 
bellows  sealed  valves  (not  shown)  will  be  employed  as  an  insurance  measure 
against  the  possible  closure  failure  of  the  Brooks  valves  and  to  establish 
the  initial  flow  conditions  of  all  gases  before  thet  enter  the  main  chamber. 
In  the  latter  process  the  gases  are  first  vented  to  the  mechanically  pumped 
side  of  the  reaction  chamber;  when  set  flow  rates  are  established,  the 
supply  is  immediately  switched  to  the  already  heated  cha.  Ker. 

The  chamber.  Figure  3,  is  composed  of  a  water  (or  He)  cooled,  vacuum 
tight  quartz  chamber  mounted  on  a  stainless  steel  plate.  Gas  entering  the 
chamber  will  be  deflected  to  the  hemispherical  cap  to  produce  a  more  evenly 
distributed  flow,  additional  mixing  and  additional  path  length  for  the 
establishment  of  equilibrium  flow  patterns.  The  r-f  heated,  barrel  type 
susceptor  will  be  a  nominally  hexagonal,  high  purity  SiC-coated  carbon 
device  supported  by  a  high  purity  Poco  graphite  assembly  and  a  stainless 
steel  rod  attached  to  a  hydraulic  and  motorized  assembly.  Figure  4,  which 
will  provide  the  capability  of  rotation  arui  vertical  movement  of  the 
susceptor.  The  latter  motion  will  allow  a  change  of  samples  in  a  loading 
chamber.  Figure  4,  without  exposing  the  fused  quartz  chamber  to  air. 

Pressure  in  the  chamber  will  be  measured  by  a  capacitance  manometer  which 
has  a  feedback  and  control  loop  through  a  second  MKS  controller  to  an 
automatic  butterfly  valve,  Figure  5;  thus,  control  of  flow  and  pressure 
will  be  independent  of  each  other.  Temperature  will  be  measured  initially 
by  an  optical  prometer;  however,  other  methods  are  being  exploded. 

As  indicated  above,  the  initial  substrates  will  be  high  purity  float  zone, 
n-type  Si  single  crystal  wafers  having  resistivities  in  the  range  of  240- 
400  ft-cm.  Both  (100)  and  (111)  wafers,  previously  etched  and  polished, 
in  order  to  compare  the  effect  of  substrate  orientation  on  the  efficacy  of 
deposition  will  be  employed. 


As  shown  in  Figures  3,  4  and  5,  a  set  of  three  gate  valves  will  be 
used  to  isolate  the  various  parts  of  the  reaction  and  pumping  sections. 
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The  di f fusion- pumped  side  will  be  used  to  evacuate  the  total  system  and 
to  differentially  pump  an  RGA  (not  shown).  The  other  pumping  assembly 
includes  two  cold  traps  to  freeze  HC1  before  it  reaches  the  automatic 
valve  and  mechanical  pump.  (There  is  some  concern  that  SiH^  may  also  be 
frozen  out  to  dangerous  concentrations;  thus,  this  aspec  ■  of  the  procedure 
will  be  closely  monitored.)  The  HC1  basin  has  been  replaced  by  a  particle 
trap  located  above  the  first  cold  trap.  The  upper  manual  butterfly  valve 
allows  a  larger  opening  in  the  automatic  valve  and  therefore  closer  control 
of  the  pressure.  A  molecular  sieve  is  also  included  to  prevent  mechanical 
pump  oil  from  backstreaming  into  the  valves  and  traps. 

In  conclusion,  this  system  has  evolved  around  the  investigators'  needs 
for  very  tight  control  on  flow  rates,  pressure,  temperature  and  purity. 
Numerous  people  at  laboratories  around  the  world  have  discussed  their 
experiences  and  offered  suggestions  on  CVD  design;  many  of  these  ideas  are 
included  in  these  drawings.  It  will  continue  to  evolve  until  the  objective 
of  the  reproducible  production  of  single  crystal  6-Si C  thin  films  is 
attained. 

III.  R.  F.  Sputtering 

An  alternative  method  of  depositing  SiC  is  via  r-f  glow  discharge 
sputtering.  In  this  process  the  target  plate  of  material  to  be  deposited 
directly  or  reacted  with  a  gas  to  produce  a  synthesized  film  is  connected 
to  a  negative  r-f  voltage  supply.  The  substrate  and  holder  face  the  target 
and  the  former  may  be  biased,  heated,  cooled  or  some  combination  of  these. 

A  gas,  usually  Ar,  is  subsequently  introduced  as  a  medium  in  which  a  glow 
discharge  can  be  initiated  and  maintained.  As  a  result  of  the  plasma, 
positive  ions  strike  the  target  plate  and  remove  principally  neutral  target 
atoms  by  momentum  transfer,  and  these  condense  into  thin  films. 

Two  approaches  are  planned  at  this  point:  1)  reactive  sputtering 
wherein  a  high  purity  Si  substrate  is  sputtered  in  a  reactive  gas  of 
methane  and  Si  and  C  species  combine  on  a  heated  substrate  to  form  SiC 
and  2)  direct  diode  sputtering  of  a  high  purity  SiC  target  onto  a  cooled 
or  heated  substrate.  It  should  be  noted  that  in  the  first  approach  the 
Si  target  may  become  conditioned  with  C  such  that  a  thin  layer  of  SiC  is 
produced  and  the  process  approach  that  of  the  second  method. 


^•*,T*r*  ■■ 
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Reactive  sputtering  or  sputtering  of  a  target  of  the  compo¬ 

sition  of  the  desired  epitaxially  grown  crystal  have  produced  reasonably 
good  films  of  g-SiC.  Wasa  et  al,  ^  have  prepared  both  amorphous  and 
polycrystalline  u-SiC  films  by  sputtering  one  inch  pressed  disks  of  .1 
and  g-SiC  of  approximately  99.7“  purity  onto  substrates  i  f  glass, 
sapphire  and  (111)  silicon.  The  polycrystalline  films  were  only  formed 
when  the  substrates  (except  glass)  were  heated  to  >  700’C.  The  most 
interesting  aspect  of  this  research  was  that  the  deposited  u-SiC  trans¬ 
formed  to  g-SiC  during  annealing  in  vacuum  at  1100°C  for  one  hour. 

1 8 

Berman  et  al_,  have  also  reported  a  method  of  converting  amorphous  Si C 

into  crystalline  SiC  by  annealing  at  high  temperature.  They  also  used 

pressed  disks  of  a-SiC  as  the  target  for  sputtering.  The  reported 

electron  diffraction  patterns  taken  in  reflection  appeared  to  include 

19 

polycrystalline  rings.  Learn  and  Haq  investigated  the  preparation  of 
single-crystalline  films  of  g-SiC  on  a-SiC  substrates  by  reactive  sputter¬ 
ing  in  acetylene  gas  using  elemental  Si  as  the  target  material.  Finally, 

20 

Nishino  et  aK  have  produced  their  own  B-SiC  target  from  crystals  etched 
from  a  silicon  solution  and  ground  into  powder.  A  (111)  silicon  wafer  was 
again  used  as  the  substrate  and  heated  to  1000°C  and  above.  These  investi¬ 
gators  concluded  that  the  substrate  temperature  should  be  kept  at  around 
1200°C  to  obtain  good  single  crystal  coatings.  The  sputtered  films  were 
oriented  in  the  same  [111]  direction  as  the  silicon  face. 

In  our  own  work,  we  are  currently  modifying  rather  extensively  and 
upgrading  an  older  sputtering  unit.  Feedback  flow  and  pressure  controls 
similar  to  that  employed  in  the  LPCVD  work  will  also  be  used  in  the 
re modeled  sputtering  system.  We  are  currently  awaiting  arrival  of  spare 
parts,  new  equipment  and  targets. 
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APPENDIX  I.  SOLGASMIX-PV  LISTING 


•  •rm.i.t,  r,c. 

c 

c  •••  solcassix-pp 

c  _ 

c 

C  CALCULATES  EQOILlrSIA  AT  COBSTAMT  PRESSURE  OB  9OL0RE.  THE  PBOC.RAB  TS  A 
C  flODIflCATIOK  Of  SOLCASB21  (G.  ERIPSS08,  CHEHICA  SCBIPTA.,  8  (1975)  100). 


1 BPL1CIT  REAl*8  (A-B.O-Z) 

COHH08  4(99,10),  A  0  ( 99 , 1  0)  ,  AKT(99>,  AKTP  (99)  ,  B(10.«>9),  B0(10). 
SG  (99) ,  KH  (20) ,  L,  B,  81,  BA,  R8,  RP  (20) ,  BL(20),  BP,  BS,  H0(99), 
S  VP,  (1  (20) ,  trot.  T  ,  TEXT (99) ,  1  (99),  TT(99),  TT  OT (20),T,NV,NVP 
DI HE K5I08  rer(99),  CE(99,6),  HE(99),  IEl(IO),  IPH(IS),  1GT  (6) , 
*108 (99),  XI (99) ,TIT1E(10) ,EL  (10) , G I  (99, 6 ) 

H  «  CLOG (10.) 

18=5 

10UT*6 

nrd)  -  1 

26  ptot  -  i. 

DO  4  8  8  *  1,  20 
98  KH  (K )  «  0 

DEAD  (18,240)  (TITLE  (I)  ,1*1,  10) 

240  POBBAT  (10A8) 

BHITE  (IOOT.2JO)  (TITLE  (I),  1*1,10) 

230  POHBAT  (• 1 • ,20X, 10AB) 

BEAD  (18,102)  L,  RP,  RB,  (RL (R) ,  8*1,  HP) 

102  PORRAT  (4012) 

BEIT  E  (I OUT, 250)  L,RP, RB, (HL(B) ,H=1,BP) 

250  EORRAT  ('080.  OP  ELERENTS  «  ',I2,5X,*80.  Of  HIXTORES  ■  *,X2.5X, 

1 ’ BO.  OP  IBVABIABT  SOLIDS  *  •,I2/,,0BO.  OP  SPECIES  PER  HIXTURE  * 

2 10  (12,  *,  ') ) 

xr  (RP  .EQ.  1)  CO  TO  146 
DO  132  8  *  2,  HP 
HP  (8)  *  HL(H-I)  «  1 
132  RL (8)  •  BL(H-I)  ♦  BL  (H) 

146  Hi  -  BL (HP)  ♦  1 
HS  *  RL(HP)  ♦  RB 
BEAD  (XR, 260)  (El  (I) ,1* 1 ,1) 

260  rORHAT  (10A8) 

DO  33  X  -  1,  as 

XI (I)  »  0. 

BEAD  (18,280)  (A  (I,J),J=1,L) 

280  POBBAT  (16F5.0) 

DO  33  J  *  1,  L 
33  A0(I,J)  -  A  (I , J) 

BEAD  (28,200)  89, T 
200  I'ORR  AT  (II, 9X, El  0.0) 

BEAD  (18,102)  (1  EL  ( J) ,  J  *  1,  L) 

BEAD  (IB, 102)  HOT,  BOB 
IP  (BCT  .ST.  0)  CO  TO  14 
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ir  (HO*  .Gl.  0)  HEAD  (IP, 102)  (IOX(l.),  N  •  1,  HOE) 

oo  2? 5  1*1, ns 

255  BIAU  (IN, 270)  TEXT  ( I) , H T  |I ) 

270  rOBHAT  (AB.2X,7E10. 3) 

CO  TC  J10 

1H  B  l  AD  (IP,  102)  (ICT(P).  P  «  1,  HO T) 

85  DO  265  1*1, «S 

265  BEAD  ( IH,  2  70)  TEXT  ( I ) •  (Ct (1 .  N)  , P* 1 , HCT) 

DO  266  1*1  ,n« 

DO  2  67  H  *  1 , 6 
267  CT  (1  ,  N)  *0. 

DO  266  N*1 , net 

266  GY(X,JCT(N))=GE(I,P) 

WHITE  (IOUT, 300)  (  T  E  XT  ( I  )  ,  (CT(1,P)  ,  P  =  1 , 6  ),1=1,RS) 

300  roan  at  cospecxes'  ,6x,  •  a*,  hi,  1  o>  ,1  ux,  »c ,  ipi.'d*  ,  mx.  •  r*  ,  1*1,  *p* , 
1/.99  ('O', 18. 31.  6 (E12.5, 3 X| ,/) ) 

310  WHITE  (IODT.230)  (TITLE(l)  ,1*1,10) 

WHITE  (IO0T.340)  (EL  (I)  ,  I-  1 ,  L) 

340  POBRAT  ('O'.ZOX.'SUESCHIPTS  OH  ELEHENT  AL  STRBOLS  OP  EACH  SPECIES', 
1/, *  0  *  » 18X,  10 (A8, 2  X)  ) 

DO  350  1=1, ns 

350  VBITE  (IOUT, 360)  T E IT (I) , ( A ( 1 , J) ,  J-1,L) 

360  rOBHAT  ('O', A8, 81, 10  (f5.  2,  H) ) 

42  BEAD  (IP, 103)  T 
WRITE  (LOOT ,122)  T 
122  P08H1T  IUHOT  *,  PB.2,  2H  */) 

IP  (nCT  .SO.  0)  GO  TO  82 
DO  84  1  *  1,  ns 
84  C  (I)  *  0. 

do  81  p  *  i,  net 

IP  (ICT(P)  .LT.  6)  TP  *  T**(1GT{«)  -  3)/8. 31433 
IP  (IGT(H)  .EQ-  6)  TP  *  DLOC  (TJ/8.  31431 
DO  81  I  *  1,  ns 

81  C(I)  *  C(I)  ♦  CC  (I ,  BJ  *TP 
CO  TO  8 

82  DO  49  I  *  1,  HS 
READ  (IB, 103)  FEP(l) 

49  G(I)  *  (PEP (1)  ♦  1000.-HF(I)/T)/8.  JU'33 
IP  (HO*  .EC.  0)  CO  TO  4 
DO  23  P  *  1,  no* 

I  *  20* (V) 

c(i)  -  -h*pep(I) 

DO  23  J  «  1,  L 
K  ■  IEL (J) 

23  C(I)  -  C(t)  ♦  A  ( I  ,  J )  /A  (K,J)  *  C  ( K ) 

4  WHITE  (1001,230)  (TITLE(I) ,1*1,10) 

•  BITE  (IOUT,  106)  (HIT  (I),  HP(I),  fEP(I|,  6(1),  1-1,  IIS) 

106  POUR  IT  ( 151,  5HHP298,  8X,  3HPEP,  91,  4HG/BT/C  •  ,  A8,  P12.  0. 2P12.  I)  ) 
CO  TO  8 

37  BtID  (IP, 103)  PTOT 
103  POBRAT  (8E10.0) 

8  ICAO  (IP, 102)  KTAL1 

CO  TO  (40,40,40,7,  34.42,37,85,26,1),  KV811 
90  BEAD  (IP, 102)  BPKT 

BEAD  (IN. 102)  (KH  (J) ,  J  *  1,  t) 

7  DO  164  J  -  1,  L 
■  ■  EH  (J) 
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60  TO  1 11.1 ,162,  H),  161,  HI,  163)  .  N 
HI  HUP  (IN, 10))  (b(.1,*|,  H  »  1 ,  H  PUT) 

CO  TO  HU 

1(2  HMD  (  '.N  ,  1  03)  B(J,  1) 

J  T  ( HI  KT  .IQ.  1)  CO  TO  !6« 

DO  160  H  *  2,  H PR T 
HO  B  <J,  N)  -  H  (J  ,  1) 

CO  TC  HU 

163  HEAD  (I  H  ,<  0  J  )  B(J.l),  STEP 
DO  165  H  •  2,  HPFT 

165  B  (J  ,  N)  -  H  (J  ,  N-  1 )  ♦  STEP 
H#  CONTINUE 

GOTO  (35,  39.34,}#),  mil 
.  35  DO  370  1*1  ns 
370  HEAD  (IN, 103)  1(1) 

GO  TO  166 
39  BEST  *  0. 

DO  169  J  *  1,  L 

IP  ( B  ( J,  1 )  -GT.  BEST  .AND.  KK  (J)  .IT.  4)  BEST  *  B(J,1) 

169  CONTINOE 

BEST  *  BEST/CPEOAT (ES) 
do  167  i  *  i,  ns 
167  »(I)  *  BEST 

166  do  135  n  =  1,  HP 
hot  <n)  *  0. 

ha  *  nr(R) 

HB  -  BL(H) 

DO  135  1  «  BA,  BB 
135  TTOT(fl)  «  TTCT(H)  4  J(I) 

3#  NP  *  0 
27  NP  *  NP  «  1 

DO  75  J  -  1,  ns 
75  NO  (1  )  *  1 

DO  22  J  «  1,1 

BO ( J  )  =  B  ( J, NP) 

J  -  ] EL ( J) 

If  (HH(J)  .11.  U)  GO  TO  170 
PI(J)  *  (G(I)  ♦  H*  BO  (J)  )/A  ( I  y  J) 

CO  TO  2? 

170  XX  ( T  1  *=  i<0  (J ) /  A  (1,0) 

ir  ( L 0  ( J )  .BE.  0.)  CO  TO  22 
DO  3  0  1  -  1,  HS 
IP  (A(I.J)  . Eg.  0.)  GO  TC  30 
NO (I)  *  0 
1(1)  -  0. 

30  CONTINUE 
22  CONTINUE 
C AIL  GASOl 

BUTE  (LOOT, 230)  (TITLE  (I).  1=1,10) 

■  Bill  (7 OUT,  111)  T,  PTOT 

111  EOBHAT  (UHOT  «,  P8.2,  2H  K/  4H  P  1PE10.3,  #H  ATS/) 

I r  (HV.EO.  0)  GO  TO  210 
UBITE  (I  OUT ,  220)  T 
220  FOBBA7  (#H  1  »,  1PE10. 3,  2H  1/) 

210  CONTINOE 

IP  (P  .IE.  BP)  GO  TO  6# 

BUTE  (IOOT,  116) 
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118  FORI)  AT  (/SON  THE  CQUILI  Hkl  tin  COHPOSITION  HAS  ROT  DEER  OBTAINED) 

CO  TC  28 

6«  If  ( ft  .  EQ.  0)  RRITE  (IOOT.121) 

121  rOINAT  (JJH  THE  SHALL  1-VALUES  ARE  NOT  EXACT/) 

00  17  J  -  1.  L 
IF  (RH(J)  .LI.  4)  CO  TO  17 
K  -  IEL(J) 

XI (*)  -  0. 

do  so  i  «  i,  as 

SO  XI  (K)  ■  XI  (K)  ♦  A0(I,J)»T(I) 

17  CONTINUE 

WHITE  (I  OUT, 105) 

105  FOBHAT  ( 15X,7HX*/flOLE,8X.6HT/noLE,10X,SHP/ATH,8X,8RACTt7ITY) 

DO  1*0  8  «  1,  HP 
IP  (II  .CT.  1)  WRITE  (IOUT,  125) 

125  iPORH  AT  (4  1  X,  1IHH0LE  FRACTION) 

8A  ■  HF (8) 

HB  -  HL  (H) 

DO  100  I  =  flA,  ns 

1*0  WHITE  (I0UT.12U)  TE  XT  (I )  ,  XI  (I),  1(1),  TF  (I)  ,  ART  (I) 

12*  FORflAT  (•  a,A8,4(3X,E12.5)  ) 

If  (NR  .CT.  0)  WHITE  (IOUT,(20)  (TEXT(I),  XI  (I) ,  T(I),  I  «  HI,  MS) 
120  FOBHAT  (/(«  •  ,  AS,  2  ( 3  X,  E 1  2.  5)  ) ) 

CALL  SPEQUA 

28  IP  (BP  -  R  PUT)  27,8,8 
1  STOP 
EHD 
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SUBBOUTJNt  C  A  SOL 
1 HPL1C1T  nil'B  (A-H.D-Z) 

COHBCW  4(99.10),  40(99,10),  AAT{99),  A*TP(99),  B{10,99),  bt)(10) 
SC (99) ,  KH{20),  L,  fl,  Hi,  BA,  HB,  Br  (20) ,  BL(20),  BP,  BS.  *0(99) 
S  HP,  P]  (20)  ,  PTOT,  T,  TEXT(99),  1  (99),  IP(99),  TTOT(20)  ,1,1V,  *H 
Cl  BE  kSION  r(99),  I P4S ( 20 )  ,  1SOL(20),  LI (99) ,  BSUB  (500) ,  OPI(IO) 
SB  (20 ,21)  ,  TPTOT(20),  11(99) 

BI  *  0 
171  IS  -  -1 
BSUB  *  0 
RG  *  0 

ss  ison  ■  o 

BS4  «  0 

IP  (BS  . LT.  fl 1)  CO  TO  47 
bo  52  1  =  Hi,  BS 
IP  (1(1)  . to.  0.)  CO  TO  52 
BS4  -  B54  *  1 
2SOL  (BS4)  *  I 

ISO!!  *  ISUB  ♦  2**  (I  «  HP  -  HI) 

S2  COST  1  RUE 
47  HPA  •  0 

■C  «  NC  ♦  1 

IP  (NG  .BQ.  501)  HG  =  1 
HSOH(RG)  '  I  SUP. 

TSOH  *  0. 

DO. 152  H  *  1,  BP  j  _•  ' 

IP  (ITOT(H)  -EO.  0.)  GO  TO  152 
HP4  *  BP4  ♦  1 
IP4S  (P.P4 )  *  H 

BSOH  (BG)  *  BSOfl(MG)  ♦  2*»(fl  -  1) 

TSOH  *  TSOH  ♦  YTOT(B) 

152  COHT I  NOE 

IP  (RSUB(KG)  .GE.  IS  .ABO.  BPS  ♦  BS4  .  LE.  L)  GO  TO  69 
59  IS  *=  15  ♦  2 

BO  154  I  *  1,  HPA 
H  •  1  PAS  (ii ) 

154  ITOT(H)  *  0. 

ISOH  -  0 
HPA  «  1 
TTOT(1)  =  1. 

IP  ( HSA  .EQ.  0)  GO  10  73 
DO  68  B  *  1,  HSA 
I  *  I SOL  (N) 

68  1(1)  *  0. 

HSA  *=  0 

73  IP  (IS  .EO.  1)  GO  TO  74 
IT  *  IS 
61  H  *  0 
71  H  *  H  4  1 

IP  (2**H  .IE.  IT)  GO  TO  71 
IP  (fl  .CT.  HP)  GO  TO  57 
HPA  «=  HPA  «  1 
1PAS  (HPA)  ■  H 
TTOT(H)  -  1. 

HA  *  BP  (fl) 
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TO  59 
«  DC  ♦  1 


rib  •  nL(fl) 

DO  150  I  «. HA.  HO 
ir  «  mo ( x )  . CO-  1>  t  (1)  ■  *s,,n 
150  CONTINUE 
CO  TO  136 

SI  I  •  m  «M«f)  -  up 
If  (I  .GT.  HS)  RETURN 
If  (»0(I)  .EC.  0)  GO  TO  59 
nsi  ■=  hsa  ♦  1 
ISOl  (USA)  *  I 

isun  *  isun  ♦  2»* (n  -  1) 
t(D  =  isua 

136  IT  -  IT  -  2**  (H  -  1) 

I r  (IT  .CT.  1)  GO  TO  6 1 
ir  (HPA  »  USA  .GT.  t)  GO 
7u  If AS  ( 1)  -  1 

If  (NSUB (MG)  .GI.  IS)  _ 

If  (SG.E0.501)  8G=1-S^ - 

NSWH  (NG)  *  IS 

69  If  (»G  . EO.  1)  CO  TO  129 
DO  1 48  K  -  2,  NG 

If  (NSUH(NG)  .EO.  NStm(K-H)  GO  TO  59 
148  CONTINOE 

if  (NSoa(ac-i)  .2Q.  nsoa)  nsuh(hg-I)  *  -nsuh(nc-I) 

If  (NG.LT.  3.0B-BSUn(RG-2).  NE.-NSUB.  03.NSUH(HG)  .  EG. HSUS) 
Nsun  (ng-2)  *  -nsuh(Ng-2) 

129  DO  142  B  *  1,  HP 
HA  *•  Hf  (8) 

HB  =  HL(H) 

If  (ITOT(B)  .GT.  0.)  GO  TO  130 
DO  126  I  =  BA,  HB 
AST  (1 )  »  0. 

A  RTF  (I)  =  1. 

1(1)  «  0. 

126  tr(I)  *  0. 

GO  TO  182 

130  DO  45  I  *  BA.  HB 

If  (HO(I)  .EQ.  1  .AND.  X  (I)  .LT.  1.  E-8)  I  <t)  *  1-E-8 
AKTf(I)  *  1. 

45  LX(I)  =  0 
CALL  A BE 8 

If  (ITOT(fl)  -EO.  0.)  GO  TO  47 
142  CONTINUE 

LSI  »  L  ♦  BP  A  ♦  HS  A 
LS  ®  LSI  -  1 
LS2  =  LS  ♦  2 
134  DHIN  >*  1.E-6 
I T  A  8  »  0 

I  TAR  J  *  BL  (1)  -  HS 
16  DO  6  0  -  1,  LSI 
DO  6  K  -  J,  LS2 
6  R  (J,  K)  *  0. 

DO  9  N  -  1,  BPA 


LI  «  L  ♦  N 
n  «  IfAS(N) 
HA  -  Hf (B) 
HB  -  HL  (H) 


GO  TO  129 
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do  s  j  *  AA,  AH 

IF  (T  (1)  .  il>.  0.)  CO  TO  5 

F(I)  '  C  (1  )  ♦  DIOC  (ART  (I  )  ) 

R(11,L52)  »  RfL1,ES2)  ♦  P(J)*t(I) 

DO  77  J  c  i,  l 

IF  (A(I,J)  .  EQ.  0.)  CO  TO  77 

AT  -  A  (I.J)  *t  (I) 

R  (J, LI)  *  mJ,D)  ♦  AT 
R  (J,  IS  2)  =  K(J.LS2)  «  AT*r(I) 

DO  BC  It  »  J,  L 

80  P(J,R)  *  R  ( J  ,  R)  ♦  AT  *A  ( I  ,  K) 

77  CONTINUE 
5  CONTINUE 

DO  9  J  *  1  ,  L 

9  P  ( J ,  L  5  2 )  *  H  (  1, LS2)  -  R(J,E1) 

IP  (ASA  „EC.  0)  CO  TO  63 
DO  67  N  *  1,  ASA 
I  =  ISOL(N) 

K  *  1  4  APA  «  N 
R  (A,  LS2)  =  C  (I) 

DO  67  J  »  1,  L 
67  R<J,  R)  *  A  (I.J) 

63  DO  31  J  =  2,  ESI 
N  »  J  -  1 
DO  31  It  =  1,  N 
31  R  (J ,  K)  =  R  (K,J) 

ir  (AT  . EC.  0)  GO  TO  156 
DO  131  J  ■  1,  L 
131  B0(J)  =  b(J,NP) 

DO  172  N  =  1,  APA 
A  *  JPAS(N) 

IF  ( AF  (A)  .NE.  AL(A))  GOTO  172 
AA  «  A  f  (  A ) 

DO  89  J  *  1,  E 

IF  ( A  (AA , J )  ,E0.  AO(AA.JJ)  GO  TO  89 
AKVOT  =  (1.  -  A  (AA,J)/A0  (AA,  J))  *1  (BA) 

DO  91  K  -  1,  L 

IF  (A(AA.K)  .EQ.  A0(AA,K))  BO  (R)  *  BOJK)  ♦  AKTOT*  A  (A  A ,  R) 
91  CONTINUE 
GO  TO  172 
89  CONTINUE 
172  CONTINUE 
156  DO  Oil  R  x  1,  L 

ir  (KH(K)  .IT.  4)  GO  TO  4U 
DO  83  J  =  1,  LSI 

83  R(J.ES2)  *  R(J,LS2)  -  PI(K)*i(J,E) 

«*>  R  (K#  LS2)  *  R  (K  ,  LS2  )  4  BO  (A) 

DO  10  It  «  1,  IS 

IF  ( FH  (K)  .61.  3)  GO  TO  10 

ELBAX  =  0. 

DO  11  J  *  K,  ESI 

IF  (ABS (B { J. R) )  .IE.  ELBA!  .OR.  KH(J)  -CT.  3)  CO  TO  11 
HBOV  »  J 

ELAAX  *DABS (P (J.K)  ) 

11  CONTINUE 

IP  (ELAAX  .61.  0.)  GO  TO  36 

IP  (N  .CT.  L  .OR.  B0(K)  .EQ.  0.)  60  TO  10 
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CO  TC  59 

n  ir  <nntm  . t g.  k>  cc  io  u 
DO  15  H  «  K,  LS2 
sadbit  •  R  (naow, n) 

R  (BN  CB, R)  “ 

15  i  (R,  8)  *  RAOBTT 
1J  IA  >  K  •  1 

DO  46  J  ”  KA,  LSI 
RK70T  ■  B ( J, K)/R (K ,K) 

DO  16  H  •  KA ,  LS2 
46  R(J.A)  «  R(J/M)  -  RKVOT*R(K,N) 

10  CONTINUE 

DO  20  R  *  1,  LSI 
K  -  IS2  -  M 

If  (KH  (K)  -GT.  3)  CO  TO  20 

ir  (R(K.K)  .BE.  0.  .AND.  B(K,LS2)  -  BE.  0.)  CO  TO  6. 

Pt  (K )  -  0. 

K  •  K  -  L  -  nPA 
I r  (K  -LE.  0)  GO ' TO  20 
I  «  ISOL  (l. ) 

I  (I)  =  0. 

CO  TO  55 

62  VI  (K )  “  B (K,IS2)/R  (K.R) 

KA  »  K  -  1 

If  (KA  . EQ.  0)  GO  TO  20 
DO  5t  J  <=  1,  KA 

58  B(J,LS2)  =>  8  (J  ,  LS2  )  -  PI(K)«8(J.K) 

20  CONTINUE  „  ..  „„  .. 

If  ( I  VAR  -BQ-  0  .OB.  IVARJ  . GE.  0  .OR.  SLAB  .LT.  0.1)  GO  TO  66 

DO  7  C  J  *  1 ,  L  _ 

IT  (DABS  (PI  (J) ) -GT.  1.E— 8.AND.DABS(OPI(J)/PI(J)-1-).GT. DBIN) 

ICO  TC  65 
70  CONTIHUE 

aa  »  o 

ir  (NG  .EQ.  1)  GO  TO  155 
DO  157  K  =  2,  NG 

IT  (SSOB(NG)  -EQ.  -BSUB(K-I))  NB  *  HB  ♦  1 
157  CONTINUE 
155  DO  143  R  «  1,  BP 

If  (TTOT(B)  .GT.  0.)  GO  TO  143 
CALL  IBEN 
TfTOT(B)  -  0. 

DO  144  I  =  BA,  B0 
TfTOT(B)  -  IFTOT(B)  *  T  f  (I) 
t*(I)  =  Tf(I) 

ART (I)  *  0. 

144  If(I)  =  0. 

IP  (NV.EQ.O)  GO  TO  3  26 
T  A*  0  . 

DO  325  I*  1  ,6L  (1) 

325  TA  *T  (I) *  TA 

PTOT  *T  A*.  0  821 *T/7 

326  COBTXHUE 

IT  (B  . EQ.  1)  TfTOT(B)  =  I f TOT  ( B) /PTOT 
143  COB? INOE 

145  oirn  *  1. 

DO  158  R  «  1,  BP 
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IP  ( 1T0T  (H )  .CT.  0.  .08.  IPTOT(fl)  . LE.  DIFB)  CO  TO 

k  a  «  b 

dipb  *  y r t ot <n) 

1SB  CONTINUE 

ip  (difb  .  f q.  i. j  co  to  i 3 o 

IT  (88  .Ey.  Cl)  GO  10  159 
NS  *  S8  -  1 
TPTOT(KA)  =»  1 . 

CO  TO  1  <1  S 

159  bsub  «  nsdbjng) 

TTOT(KA)  =  1. 

BA  *  BP(KA) 

BB  ®  BL(KA) 

DO  153  1  -  BA,  H  8 
153  1(1)  =  Y  S  Ofl*  X  I  (I )  /  Y  PTOT  (KA) 

CO  TC  87 

138  If  (PIS  .  LT.  81)  CO  TO  51 
DIPB  «  0. 

DO  58  I  «  81,  BS 

IP  (80(1)  .tQ.  0  .OB.  T  ( I)  .GT.  0.)  GO  TO  58 
TIB  *  -G(l) 

DO  56  J  =  1,  L 
56  PI  A  *  PI  A  ♦  A  (I  ,  J)  *P1{J) 

I?  ( PI A  .IT.  DIP B)  CO  TO  58 
KA  c  I 

Dim  »  pi  a 

58  CONTINUE 

ir  (U1PB  . EQ.  0.)  GO  TO  51 
1  (Kt )  «=  TSUB 
CO  TO  55 

51  IP  (BX  . Ey .  1)  GO  TO  93 
DO  168  N  *  1,  BPA 
B  =  IPAS (N) 

IP  ( BF  (B )  .  N £.  BL  ( H)  )  GO  TO  168 

BA  *  BP ( B) 

DO  173  J  =  1,  L 

IF  ( A  (BA, J )  .EO.  A 0  (BA , J) )  GO  TO  173 
BX  *  1 
GO  TO  171 
173  CONTINUE 
168  CONTINUE 
93  ITARJ  =  0 
GO  TO  66 

65  DO  (0  J  >  1,  1 
60  OPI (J)  *  Pi (j) 

66  SLAB  *  1. 

DO  12  N  a  1,  BPA 
LI  *  t  ♦  N 
B  «=  IPAS  (N) 

KA  *  BP  (B) 

BB  *  BL ( B) 

DO  12  I  *  BA,  BB 

IP  (T  (I)  . EQ.  0.)  CO  TO  12 

PI A  *  P(I)  -  PI  (LI ) 

DO  11  J  =  1,  j. 

19  PI  A  *  PI  A  -  A  (I  ,  J)  *PI  (J) 

11(1)  =  PI  A*  I  (I) 
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lr  (PI*  _GT.  SLAN.)  SLAfl  ■  PI  * 

12  CONTINUE 

!r.  r,tA'’  *  0-W(SlA«l  -  0.i»/(SL*N*SL*n) 

IP  (F.SA  .  EQ.  0)  CO  TO  72 
DO  S3  N  ■  1,  ns* 

I  *  ISOL(N) 

*  ■  l  ♦  bpa  ♦  h 

IP  (PI(K)  .  LT.  0.)  LI  -  0 

IP  ( I  V*B  .EQ.  0  .03.  P I  ( A)  .GT.  - T  (I)  )  CO  TO  5) 

IP  (SLAfl  -LT-  1.  .AND.  -PI(K)/TSUN  .LT.  1.E6)  CO  TO  53 
I <IJ  ■»  0. 

GO  TC  55 

53  I  (IJ  rDABS  (PI  (K)  ) 

72  TSUH  x  0. 

DO  1^8  a  -  1,  HP* 
n  «  1 FAS ( N ) 
ba  *  BP(n) 
bb  =  ni(nj 
do  29  i  a  n*.  ns 


»  0. 


47 


-KO-  50)  DUIN  *  100.  •DHIN 

l7ARJ,  =  i<ITiBj\01*OB*  11  *B°*  °  *°R*  St*n  *1‘T*  1*)  00  T0  1fc 
IF  (IVARJ  .BQ.  10)  GO  TO  88 
DO  3  N  *  1,  RPA 
B  -  IF  AS  (  N) 

ba  *  BP(.n) 

BB  =  HL (H) 

DO  3  I  «  n*,  (IB 

3  CONTINUE  °*  **HD-D»BS<rX(I))/T{I>  *01.  1.E-6)  CO  TO  16 

88  DO  149  N  =  1,  nr* 

B  *  IFAS(H) 

CALL  XBER 

IP  (BA  -BE.  Rfl)  GO  TO  1J3 
DO  92  J  =  1,  L 

92  C0NTINUe'J>  *NE*  A0<"A'J5,  l<"*)  "  MBA,J)/AO(R*,J)*r(BA) 

133  DO  149  I  x,  HA,  n B 

IP  (*0(I)  .BC.  0  .08.  T  ( I)  .GT.  0.)  CO  TO  149 
Id)  =  TTOI(B)  •!?(!) 

f!'111’  -lI<  1.E-20  .08.  LI  (I)  .  EQ-  1)  CO  TO  149 
LI(I)  X  1 

GO  TO  13i 
149  C0NTI8UE 

IP  (IVA8J  .EQ.  10)  B  -  0 

BETUBB 

END 


IP  (1(1)  -EQ.  0.)  GO  TO  29 
HI)  *  1(1)  -  SLAn*TX(I) 

IP  (1(1)  .LT.  1.E-20)  III) 
29  CONTINUE 
CALL  ABER 

IP  (ITOT(B)  .EQ.  o.)  GO  TO 
128  ISOB  =  TSUB  ♦  TTOT(fl) 

17  AR  x  X  V A  R  ♦  1 

IP  ( I VAR  .EQ.  25  .OB.  ITAH 

IF  ( I V AR  . EQ.  75)  CO  TO  59 


i 

I 


i 


I 

I 

I 

! 

I 


SUBROUTINE  A Blf 
IBPLJCIT  REAL’S  (A-H.0-2) 

COBBCH  A  (99,  10),  A0(99,10),  ART(99),  AKTP  (99) ,  6(10,99),  BO(IO) 
SG  (99) ,  KH(20),  t,  fl,  N1.  BA,  JIB,  Hr  (20)  ,  flt  (20)  ,  BP,  BS,  N0(99) 
SNP,  PI  (20)  ,  PTOT,  T,  TEYTJ99).  T  (99)  ,  TP  (99)  ,  HOT  (20)  ,V  ,  N  V,  NVV 
ITOT  |B)  •  0. 

DO  2  I  *  BA,  BB 
2  ITOT (B)  =  I1C1  (B)  ♦  Y (1 ) 

IF  (TTOT(B)  . GF..  1.1-6)  GO  TO  151 
TTOT(fl)  *  0. 

RETURN 

151  CONTINUE 

IF  (Mf.EO.  0)  GO  TO  3 26 
YA«0. 

DO  325  1  =  1. Bill) 

325  I A  =  T  (I )  «  Y  A 
PTOT=TA*.0821*T/V 

326  CONTINUE 

IP  (B  -EQ.  1)  TTOT(I)  *  TTOT  (1)  /PTOT 
DO  127  I  =  BA,  BB 
127  TP  (I)  ’  T ( I) /TTOT (fl) 

CALL  I  ACTON 
DO  191  T  *  BA,  BB 
141  ART  (I )  *  A  KTE  (I )  •  TP  (I) 

RETURN 

END 
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SUBBOUTI HE  X  BEB 
IHPLICJT  REAL'fl  (l-H.O-Z) 

C CRfl OH  *(99,10),  * 0(99,10),  *KT(99),  *KTP|99),  BJ10.99),  B0(10t, 
tC  (99 )  ,  KH  (20)  ,  L.  A,  n  1 ,  HA.  R  B ,  HP  (20)  ,  Hi  (  20)  ,  HP,  I1S.  H  0  (99)  , 
*«P,  PI  (20) ,  PTOT,  T,  TEXT  (99) ,  X (99) .  IP(99],  TTOT (20) , X, BV, HVV 
D1 A2 1SIOB  01  P  (99) 

PLOG  -  DHAI1  (DLOC(PTOT)  ,  10.) 

H*  -  HP  (H) 

HB  «  BL ( H) 

00  38  I  *  HA,  HB 

XP  ( HO  (I)  .KQ.  0  .0*.  T  ( I)  -fit.  0.)  GO  TO  38 
PI*  »  -0(1) 

DO  87  J  =  1,  L 
87>  PI*  *  Pt*  ♦  A (I ,  J) •PI  (J) 

IP  (PI*  -CT.  FLOG)  PI*  =  PLOG 
*KT ( I)  -DEXP  (PI*) 

TF  (1 )  =■  AKT  (I) 

38  COHTJHOe 
I VAC  >  0 

1«7  IVA8  *  IX*  R  ♦  1 

IP  (IXAR  .  LT.  25)  CALL  P*CTOB 
00  139  I  m  HI,  HB 
OIF(l)  *  T  P  (1) 

139  TP  (I  )  =  AKT(I)/AKTJ  (I) 

DO  137  I  =  HA,  RB 

IP  (TP(I).GT.O..AND.D*BS{OTF(I)/TP{I)-1.).GT.1.E-«)GO  TO  H7 
137  COHTIHOE 
BETH  BX 
END 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  FACTOR 
1 BPL1C1T  NEAL'fl  (A-H.O-t) 

CCnnCH  A  (99,  10)  ,  *0199,10).  ART(99),  AKTF(99),  8(  10,99),  B0(10), 

Sc  (99 ) ,  Mi(20),  t,  n,  ni.  ha.  bb,  nr  (20) ,  m,(20),  ni>,  ns,  *0(99), 

SUP,  PJ  (20)  ,  HOT,  T  ,  TEXT  (99)  ,  T  (99)  ,  IP(99),  TTOT  (20)  ,  V  ,  H  V ,  NVY 

THE  ACTIVITY  COEFFICIENT  EXPRESSIONS  rOB  THE  SPECIES  IN  THE  NONIDEAL 
HIXTURES  ARE  INSERTED  A rT ER  THE  STATEnENT  NURHER  THAT  COBPESPOHDS  TO 
THE  B1 X  TUBE  NUHBER.  POH  A  SOLID  WITH  VARIABLE  STO > C HI  On ETR I  WHICH  IS 
BEGIRDED  AS  A  SEPARATE  niXTURE,  THE  COHPOSITION  HOST  ALSO  BE  SPECIFIED 
(RAXlnUn  ONE  STOICHIonETRIC  CO  EE  FI C IERT  IS  ALLOWED  TO  DEVIATE). 

NOTICE  THAT  AN  AEOUNT  OP  A  SPECIES  WHICH  IS  LESS  THAN  1.E-2D  HOLES 
IS  SET  EQUAL  TO  ZERO. 

AO  =  A  DEVIATING  STO ICH I  0 HETR I C  COEFFICIENT 

AKT  -  ACTIVITT 

AKTF  =  ACTIVITY  COEFFICIENT 

T  *  AHOUNT  OF  SUBSTANCE 

TP  =  HOLE  FRACTICN  OR  PARTIAL  PRESSURE 

GO  TC  (1,2, 3,4, 5, 6, 7, 8, 9,10) ,  ft 

1  CONTINUE 
RETURN 

2  CONTINUE 
RETURN 

3  CONTINUE 
RETURN 

4  CONTINUE 
RETURN 

5  CONTINUE 
RETURN 

6  CONTINUE 
RETURN 

7  CONTINUE 
RETURN 

8  CONTINUE 
RETURN 

9  CONTINUE 
RETURN 

10  CONTINUE 
RETURN 
END 
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SUUHCUT1NK  SPKQUA 
IMPLICIT  «EAL»8  (A-H.O-Z) 

COMMON  *(99.  10).  *0  (99,  10),  *KT  (99)  ,  A*Tr(99).  0(10,99),  B0(10), 
SC  (99 ) ,  *H(20).  L,  n,  a  1 ,  a*.  MO.  «f(20|  ,  flL (20)  ,  «P.  AS,  NO (99) , 
SBP,  PI  (20)  ,  FTOT,  T.  TEST  (99)  ,  t  (99)  ,  ir(99),  ITOT  (20)  ,  V  ,  ti  V ,  NV  V 

S  PEQU  A  IS  A  SUBRCUT1 NE  fOH  CALCULATION  OP  QUANTITIES  NHICH  ABE 
DIBIT  ABLE  PBOfl  TNI  EQUILIBRIUM  COMPOSITION. 

RETURN 

END 
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